Optimization and performance of an optical cardio-magnetometer 
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Cardiomagnetonietry is a growing field of noninvasive medical diagnostics that 
has triggered a need for affordable high-sensitivity magnetometers. Optical 
pumping magnetometers are promising candidates satisfying that need since 
it was demonstrated that thy can map the heart magnetic field. For the op- 
timization of such devices theoretical limits on the performance as well as an 
experimental approach is presented. The promising result is a intrinsic mag- 
netometric sensitivity of 63 fT/vHz a measurement bandwidth of 140 Hz and 
a spatial resolution of 28 mm. 
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1. Introduction 

Bioniagnetometry is a rapidly growing field of noninva- 
sive medical diagnostics^. In particular, the magnetic 
fields generated by the human heart and brain carry valu- 
able information about the underlying electrophysiolog- 
ical processes^. Since the 1970s superconducting quan- 
tum interference devices (SQUIDs) have been used to de- 
tect these generally very weak biomagnetic fields-. The 
magnetic field of the human heart is the strongest bio- 
magnetic signal, with a peak amplitude of 100 pT, but 
since this is still orders of magnitude weaker than typi- 
cal stray field interference the measurement of such sig- 
nals could initially only be performed inside expensive 
magnetically-shielded rooms (MSR). Progress in medi- 
cal research in the past decade has motivated a need for 
more affordable cardiomagnetic sensors. Recently, mul- 
tichannel SQUIDs were developed that no longer require 
shielding due to the use of gradiometric configurations. 
Such devices are commercially available but are still quite 
expensive in both capital and operational costs. 

Optical pumping magnetometers (0PM) have been 
widely known since the 1960a^, and offer both high sen- 
sitivity and reliable operation for research^ and applica- 
tions like geomagnetometryS.. Since OPMs usually work 
with a near room-temperature thermal alkali metal va- 
por, they avoid the need for the cryogenic cooling that 
makes SQUIDs so costly and maintenance intensive. Our 
goal was to develop an affordable, maintenance-free de- 
vice that is both sensitive and fast enough to measure 
the magnetic field of the human heart. In order to be 
competitive with the well-established SQUIDs, a cardio- 
magnetic sensor has to offer a magnetic field sensitivity 
of at least 1 pT with a bandwidth of about 100 Hz. Fur- 
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thermore, the spatial resolution of the sensor has to be 
better than 4 cm, the standard separation of grid points 
during mapping. 

Since the cardiomagnetonietry community is mainly 
interested in one of the components of the magnetic field 
vector, one might think of using vector-type OPMs like 
the Hanle magnetometesi or the Faraday magnetometer, 
devices which operate in zero fields only^. However, these 
devices lose their sensitivity in the presence of even tiny 
field components in directions perpendicular to the field 
of interest. The broadening caused by such transverse 
field components must be kept well below the width of 
the magnetometer resonance2>, thus limiting those com- 
ponents to values below a few tenths of pT. Accord- 
ingly, optical vector magnetometers cannot be used for 
cardiomagnetonietry in a straightforward way since the 
heart field features time-varying transverse field compo- 
nents on the order of 100 pT. We have therefore con- 
centrated on the Mx 0PM, which exhibits a fast re- 
sponse and which has been shown to be sufficiently sen- 
sitive in an unshielded environment. Furthermore, lamp- 
pumped Mx OPMs were used for the first biomagnetic 
measurements^ with optical magnetometers in the early 
1980s, although that work was discontinued. Instead of 
lamps, we use diode lasers as a light source in order to 
build a device that will scale to the many channels needed 
for fast mapping of the cardiomagnetic field. 

2. Principle of scalar OPM operation 

Optically pumped magnetometers operate on the prin- 
ciple that the optical properties of a suitable atomic 
medium are coupled to its magnetic properties via the 
atomic spin. The ensemble average of the magnetic mo- 
ments associated with the spins can be treated as a classi- 
cal magnetization vector M ~ Npgp IJ-b^P) /^ in space. 
Here Nf{P) = NptT:{pF'P) with iv pp = 1 is the total 
angular momentum of Np atoms in an optical hyper- 
fine level F where pp is the density matrix and gp the 
Lande factor of the state. Optical magnetometers detect 




Fig. 1. Basic geometry of the Mr-magnetometer setup: 
The laser (la) emits a beam that traverses the sample 
(sa) at angle 9 with respect to the magnetic field Bq. 
The transmitted power is detected by a photodetector 
(pd). The static magnetic field Bq is aligned along the 
^-direction. The oscillating magnetic field Bi is aligned 
along the x-direction. 



changes of the medium's optical properties induced by 
the precession of M in a magnetic field Btot- The fre- 
quency of this precession, the Larmor frequency lol-, is 
proportional to the modulus of Btot '■ 



UJL = — r Btot 



7F|Bt 



(1) 



For Cs the constant of proportionality, jf , has a value of 
2tt X 3.5Hz/nT. All atomic vapor magnetometers measure 
the magnetic field via a direct or indirect measurement 
of the Larmor frequency. 

A. The Mx magnetometer 

In the case of the M^ magnetometer, a magnetic- 
resonance technique is used to measure the Larmor fre- 
quency directly, by employing two perpendicular mag- 
netic fields Bo and Bi. The static magnetic field Bq 
is aligned along the z-direction. As Fig. ^ shows, the 
k- vector of the laser beam lies in the yz-plane and is ori- 
ented at an angle 9 with respect to the z-direction. The 
magnetometer is sensitive to the modulus of Bq. The 
oscillating magnetic field Bi = xBif cos Wrf i is aligned 
along the x-direction with an amplitude much smaller 
than Bq. 

In order to introduce the basic concepts we discuss the 
simplest case of an i^ = 1/2 state. The motion of M 
under the influence of Bq and Bi is then given by the 
Bloch equations: 
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The first term describes the precession of M around the 
magnetic fields. The second term describes the longitu- 
dinal (71) and transverse (72) relaxation of M. The third 
term represents the effect of optical pumping with circu- 
larly polarized light that creates the magnetization. It 
can be treated as an additional relaxation leading to an 
equilibrium orientation aligned with the k- vector of the 
incoming light at the pumping rate Tp. Both relaxations 
add up to the effective relaxation rates ri.2 = 71.2 + Tp- 
In the case of small B,-{ amplitudes, Eq. ^ can be solved 
using the rotating-wave approximation^*^ which leads to 
a steady-state solution where M rotates around Bq at 
the driving frequency Wj-f . 

The optical property used in the M^ magnetometer is 
the optical absorption coefficient which determines the 
power, P, of the light transnntted through the medium. 
For circularly polarized light, the transmitted power is 
proportional to the projection of M on the k- vector of the 
incoming light. Therefore, the precessing magnetization 
results in a modulation of the absorption index measur- 
able as an oscillation of P. The in-phase and quadrature 
components of P with respect to the driving field can be 
obtained from Eq. 0: 

Pip(<5) = -Posin(2g) ^,^ ,^''^^^ , ^^ ,and (3) 
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Here firf — ^pBrt is the Rabi frequency and S = uj-^i — ujL 
the detuning of the oscillating field Bi from the Larmor 
frequency. The constant Pq combines all factors such 
as the initial light power, the number of atoms in the 
sample, and the cross section for light-atom interactions 
determining the absolute amplitude of the signal. The 
components can be measured using phase-sensitive de- 
tection. The signals are strongest for 9 — 45°, which was 
used in all experiments. 

Both Pip and Pqu show resonant behavior near 5 = 0. 
Pqu has an absorptive Lorentzian line shape, and Pjp has 
a dispersive Lorentzian line shape with the same half 
width expressed as 



AcjHw = V"rfr2/ri 



Tl=T2VsTl. 



(5) 



Here S ~ VP' /(riT2) is the saturation parameter of the rf 
field. Figure[2a) shows measured line shapes under con- 
ditions optimized for maximal magnetometric sensitivity 
(see Sec. ^ for details). 

Signal Pip is of particular interest because it has a dis- 
persive shape, featuring a steep linear zero-crossing at 
(5 = 0. In this region Pip can be used to measure the 
deviation of Bq from the value that corresponds to Wrf. 
The same is true for the deviation of the phase difference 
Lf between the measured oscillation and the driving field 
from -90° (see Fig.jSJb)). The phase difference tp can be 
calculated from Pip and Pqu, yielding 



tan(/3 



P 



ip 



£2 
5 



(6) 




-600 -400 



-200 200 

(5/27r (Hz) 



400 600 



Fig. 2. (a) Measured magnetic resonance lineshapes of 
the in-phase (ai) and quadrature signals (02), measured 
in a single sweep of 20 s with the cardiomagnetometer 
placed in a poorly shielded room. Magnetic 50 Hz line 
interference was suppressed using a 4*'* order lowpass fil- 
ter (time constant 10 ms). The half- width, derived from a 
fit, was Acjhw/Stt — 173Hz. (b) Magnetic resonance line 
shape of the oscillation phase measured with respect to 
the driving field i?rf • The data was obtained in real time 
using a digital lock-in amplifier. The fitted half-width is: 
Awgw/27r= 109 Hz. 



The phase signal changes from 1^9 = at low frequencies 
to iy9 = — TT at high frequencies. For practical reasons it 
is preferable to shift the phase by 90° so that it passes 
through zero in the center of the resonance (5 = 0). This 
can easily be done by shifting the reference signal by 7r/2 
using the corresponding feature of the phase detector. In 
mathematical terms that 90° shift is equivalent to the 
transformation P„„ 



-Pip and Pip ^ Pqu, yielding 



tan(/3' = - — . (7) 
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The width of the phase signal Aa;|j^ is smaller than 
Awhw because it is not affected by rf power broaden- 
ing, i.e., it is independent of firf : 



^HW 



= r2 < AWHW 



(8) 



The narrower lineshape of the phase signal is exactly 
compensated by a better S/N ratio of Pip [see Eqs. Q, 
(|10|l . and (|ll|l ] resulting in a statistically equivalent mag- 
netic field resolution for both signals. However, since the 
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Fig. 3. Nyquist plots for three different values of 5*. The 
solid line is for 5* = 1, the dashed lines are for S = 0.25 
(nearly round) and 5 = 4 (elliptical) , respectively. If cj^f 
is scanned towards increasing values the system evolves 
clockwise through the Nyquist plot. 



lineshape of the phase signal depends only on r2, it is 
easier to calibrate in absolute field units. Furthermore, 
light amplitude noise, for instance caused by fluctuating 
laser intensities, does not directly affect the phase sig- 
nal, since both Pip and Pqu scale in the same way with 
light intensity. Only the much weaker coupling via the 
light shift can cause the phase signal to reflect light am- 
plitude noise. Considering those practical advantages of 
the phase signal we concentrate in the following sections 
on the sensitivity of the phase signal to magnetic field 
changes. 

B. Nyquist plots 

The hneshapes Pip {5) , and Pqu [5) of the magnetic reso- 
nance have a major influence on the magnetometric sen- 
sitivity. The magnetic resonances in Eqs. Q and Q 
can be interpreted as a complex transfer function i(w) 
connecting the current that drives the rf-coils, U — 
Iq exp^iujt) + c.c, and the photocurrent Ip = t{u))Ii of 
the photodiode. By setting the effective transverse re- 
laxation rate r2 as the unit of frequency and using the 
normalized detuning x (= 5/T2)^ t can be written in di- 
mensionless units as: 



x/5(» + x) 
*-*°l + 5 + x2 



(9) 



A parametric plot of t{x) in the complex plane — called 
a Nyquist plot — was found to be useful for the inspec- 
tion of experimental data. In this representation t{x) 
appears as an ellipse with diameters dr and di for the 
real (in-phase) and imaginary (quadrature) components 
respectively (see Fig. IS)): 
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Fig. 4. Nyquist plot of a resonance with S = 1 when 
an interfering sine wave of amphtude r^ and phase (fi is 
added. A phase offset of A(p in the demodulation due to 
a poorly adjusted lock-in phase leads to a rotated ellipse. 



The saturation parameter of the rf transition, S, can be 
extracted from the ratio of the two diameters: 
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Figure 01 shows a Nyquist plot of a resonance for a 
situation in which an interfering sine wave is added to 
the photocurrent, leading to a shifted ellipse. The am- 
plitude Vi and the phase ipi of the interference can be 
easily extracted from the Nyquist plot. A phase shift 
in the demodulation leads to a rotated ellipse. In this 
situation the spectra of in-phase -Pip (5) and quadrature 
-Pqu(^) components as a function of rf detuning appear 
asymmetric. 

By means of Nyquist plots it is easy to distinguish 
between an asymmetry caused by improper adjustment 
of the lock-in phase and one caused by inhomogeneous 
broadening. The latter causes a deviation from the ellip- 
tical shape. One model for inhomogeneous broadening is 
to assume a gradient in the static magnetic field. Since 
we use buffer-gas cells the atoms do not move over large 
distances during their spin coherence lifetime so that in- 
homogeneous magnetic fields are not averaged out. In- 
stead, atoms at different locations in the cell see different 
magnetic fields, resulting in an inhomogeneous broaden- 
ing of the magnetic resonance line. 

Figure [S] shows calculated Nyquist plots for different 
gradients of the static field Bq. The simplest model for 
such an inhomogeneity is a constant gradient dBz/dz 
over the length Iz of the cell. This is expressed by a 
convolution of the theoretical magnetic resonance signals 
t{x) [see Eq. Q] with the normalized distribution of mag- 
netic fields g{x) which, in this case, is a constant over the 
interval 
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Since g{x) vanishes everywhere except for 



(13) 




Fig. 5. Nyquist plots for different magnetic field distri- 
butions, each scaled to fit in a square of length 1. Part 
(a) shows the deviation from circular for constant mag- 
netic field distributions. The innermost trace is for an 
unperturbed resonance. The two outer traces are cal- 
culated for field distribution widths of Xg = 40 r2 and 
400 r2, respectively. Part (b) shows the deviation from 
circular for linear field distributions. The (outer) circular 
trace is for an unperturbed resonance. The other two are 
(12) calculated for distribution widths of Xg = 5T2 and 10 r2, 
respectively. 

the convoluted resonance t' is given by 

t'{x) = I t(x~x')g{x')<lx' , (14) 

which can be evaluated analytically 
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The main effect of the constant magnetic field distribu- 
tion is to broaden the resonance, to decrease the ampli- 
tude, and to make the line shape differ from a Lorentzian. 
In the Nyquist plot this is seen by a deformation of the 
elliptical trace towards a rectangular trace as shown in 
Fig. Ela). The effect is clearly visible in Fig. EJa) for 
rather large widths of the magnetic field distribution; 
in the experiment, however, the effect can be detected 
for much smaller inhomogeneities due to the large sig- 
nal/noise ratio. 

3. Experimental Setup 

The magnetometer described here was part of the device 
used by us to measure the magnetic field of the human 
hearbiiii^. The setup was designed so that a volunteer 
could be placed under the sensor, with his heart close 
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Fig. 6. Schematic of the experimental setup. Light from 
the diode laser is delivered via an optical fiber (of) to 
the experiment in the magnetically shielded room (msr). 
The light is coUimated by a lens (le) and polarized by 
a polarizing beam splitter cube (po). The beam splitter 
(bs) reflects 50% of the beam to photodiode 1 (pdl) used 
for monitoring of the initial light power. The remaining 
beam passes a quarter-wave plate (A/4) providing circu- 
lar polarized light to the glass cell (gs) that contains the 
atomic medium. Photodiode 2 (pd2) measures the trans- 
mitted light intensity. Its signal is amplified by a current 
amplifier and fed to the lock- in amplifier (lia). The ref- 
erence output of the lia drives the radio frequency coils 
(rfc). The reference frequency of the lia is controlled by a 
sweep generator (sg). Automatic control and data aqui- 
sition is done by a PC via the GPIB bus. 



to the glass cell containing the Cs sample. For moving 
the volunteer with respect to the sensor — necessary for 
mapping the heart magnetic field — a bed on a low friction 
support was used. 

The magnetometer sensor head itself was placed in a 
room with moderate magnetic shielding. The room was 
1.7 X 2.3 X 2.5 m'^ in volume shielded by a 1 mm //-metal 
layer and an 8 mm copper-coated aluminum layer. For 
low frequencies, the shielding factor was as low as 5 to 
10, whereas 50 Hz interference was suppressed by a factor 
of 150. Inside the shielded room, surrounding the sensor 
itself, three coil pairs were placed for the three dimen- 
sional control of the magnetic field. In the z-direction 
(vertical) two round 1 m diameter coils were used. To 
make room for the patient, the spacing between the coils 
had to be 62 cm, far away from the Helmholtz optimum 
of 50 cm. The two coil pairs for the transverse magnetic 
fields {x and y directions) formed four of the faces of a 
cube 62 cm on a side. All six coils were driven inde- 
pendently by current sources so that the sum and the 
difference of the currents in each coil pair could be cho- 
sen independently. This allowed us to control not only 
the magnetic field amplitudes in all three directions, but 
also the gradients dBi/di. The field components and gra- 
dients were adjusted to produce a homogeneous field of 
5 //T in the z direction. 



An extended-cavity diode laser outside the shielded 
room was used as a light source. The laser frequency 
was actively stabilized to the F ^ 4 ^ F ^ 3 transi- 
tion of the Doppler broadened Cs Di line (894 nm) using 
DAVLL spectroscopy^'^ in an auxiliary cell. The light 
was delivered to the magnetometer sensor proper by a 
multimode fiber (800 /im core diameter). After being 
coUimated, the light was circularly polarized by a combi- 
nation of a polarizing beam-splitter and a multiple-order 
quarter-wave plate. The circularly polarized light then 
passed through a glass cell containing the Cs vapor and 
a buffer gas to prevent the atoms from being depolarized 
by wall collisions. The cell could be heated to 65° C us- 
ing hot air which flowed through silicon tubes wrapped 
around the cell holder. The light power, P, transmit- 
ted through the glass cell was detected by a photodiode 
specially selected to contain no magnetic materials. A 
current amplifier (FEMTO Messtechnik, model DLPCA- 
200) converted the photocurrent into a voltage that was 
fed to the input of the lock-in amplifier. The detection 
method resulted in a noise level 5 to 20% above the elec- 
tron shot noise in the photodiode (Fig. Cj). The digi- 
tal lock-in amplifier (Stanford Research Systems, model 
SR830) demodulated the oscillation of P with reference 
to the applied oscillating magnetic field. That field was 
generated by two extra windings on each of the B^ coils 
and was powered by the analog output of the reference 
function generator contained within the lock-in amplifier. 
The built-in function generator has the advantage that it 
delivers a very pure sine wave (phase locked to the syn- 
chronization input) and its amplitude can be controlled 
via the GPIB interface of the lock-in amplifier. 

In order to record magnetic resonance lineshapes the 
lock-in amplifier was synchronized to a reference fre- 
quency supplied by a scanning function generator. The 
data measured by the lock-in (amplitudes of the in-phase 
and quadrature signals) were transmitted in digital form 
to a PC, thus avoiding additional noise. 



4. Optimization 

Although the theory of optical magnetometry is well 
knowrA, predictions about the real performance of a 
magnetometer, especially when it is operating in weakly 
shielded environments, are difficult to make. The perfor- 
mance depends on laser power, rf power, cell size, laser 
beam profile, buffer-gas pressure, and the temperature- 
dependent density of Cs atoms. The size of the cells and 
the buffer gas pressure were dictated by the available 
cells: We used 20 mm long cells with 20 mm diameter 
including 45 mbar Ne and 8 mbar Ar with a saturated Cs 
vapor. Since the cell is oriented at 45° with respect to Bq, 
the transverse spatial resolution was 28 mm. The cross 
section of the laser beam was limited by the 8 mm aper- 
tures of the optical components (polarizers and quarter- 
wave plates). 
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Fig. 7. Root power spectrum of the photocurrent when 
the driving field is in resonance with the Larnior fre- 
quency. The data sample was recorded at 54° C under 
conditions optimized for maximum magnetometric sensi- 
tivity with a resolution bandwidth of 1 Hz (sampling time 
1 s). The amplitude measured by the lock- in amplifier 
corresponds to the upper horizontal line. The amplitude 
of the central peak is depressed, since it is slightly broad- 
ened by the Hanning window used by the FFT spectrum 
analyzer (see text). The level pi is the shot-noise level 
calculated from the DC-photocurrent. The dashed line 
marks the rms noise measured at 23 kHz. The i?sN with 
respect to the calculated shot- noise level is 5 x 10^. The 
rms noise is a factor of 1.55 higher than pj resulting in a 
i?SN of 3.2 X 10^. 



A. Intrinsic resolution 

Our magnetometer produces a signal which was propor- 
tional to the magnetic field changes. The noise of the 
signal in a perfectly stable field therefore determines the 
smallest measurable magnetic field change, called the 
noise equivalent magnetic field (NEM). The NEM is given 
by the square root, pB, of the power spectral density, p^, 
of the magnetometer signal, expressed in T/Hz^". The 
rms noise, cs , of the magnetometer in a given bandwidth 
/bw is then 

(^B ^ PB\/lh^- (16) 

A straightforward way to measure the intrinsic sensitivity 
would be to extract the noise level from a sampled mag- 
netometer time series via a Fourier transformation. How- 
ever, that process requires very good magnetic shielding 
since the measured noise is the sum of the magnetic field 
noise and the intrinsic noise of the magnetometer. Many 
studies under well-shielded conditions have been carried 
out in our laboratory, leading to the result that optical 
magnetometers are in principle sensitive enough to mea- 
sure the magnetic field of the human heart. However, 
the shielding cylinders used in these investigations were 
too small to accommodate a person. The present study 
investigates which level of performance can be obtained 
in a weakly shielded environment with a volume large 



enough to perform biomagnetic measurements on adults. 
In the walk-in shielding chamber available in our labo- 
ratory the magnetic noise level was about one order of 
magnitude larger than the strongest magnetic field gen- 
erated by the heart. In order to compensate for this the 
actual cardiomagnetic measurements were done with two 
magnetometers in a gradiometric configurationli. How- 
ever, the optimal working parameters where determined 
for a single magnetometer channel only. 

Since all time series recorded in this environment are 
dominated by magnetic field noise, the straightforward 
way of measuring the intrinsic noise could not be ap- 
plied. As an alternative approach a lower limit for the 
intrinsic noise can be calculated using information the- 
ory. The so-called Cramer-Rao lower boundsA^ gives a 
lower limit on how precisely parameters, such as phase or 
frequency, can be extracted from a signal in the presence 
of a certain noise level. For the following discussion we 
assume that the signal is a pure sine wave affected by 
white noise with a power spectral density of p^ . We de- 
fine the signal-to-noise ratio -Rsn as the rms amplitude, 
A, of the sinusoidal signal divided by the noise amplitude, 
(7, for the measurement bandwidth, /bw: 



-Rsn — — — — ■?== 

O" PV/bw 



(17) 



For a magnetometer generating a Larmor frequency 
proportional to the magnetic field, Eq. (^, the ultimate 
magnetic sensitivity is limited by the frequency measure- 
ment process. The Cramer-Rao lower bound for the vari- 
ance, K,, of the frequency measurementii is used (Ap- 
pendix ITf to calculate pb 
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For cardiac measurements a bandwidth of /bw — 100 Hz 
is required. This together with a typical value for 
i?SN of 10'* results in a magnetic field resolution of 
315 fT/Hz^/^. In order to be competitive with SQUID- 
based cardiomagnetometers that feature an intrinsic 
noise of 5. . . 20 fT/Hz ' this level of performance is not 
sufficient. For that reason we have concentrated on a dif- 
ferent mode of operation where the phase signal is mea- 
sured by digital lock-in detection. 

In this mode of operation Wrf has a fixed value near the 
Larmor frequency. The information about the magnetic 
field is obtained from the phase shift of the magnetometer 
response at that frequency. The Cramer-Rao bound for 
a phase measurement of a signal with known frequency is 
used in Appendix |H1 to calculate the NEM for that case: 
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Equations H18(l and (|19|l define the bandwidth: 
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Fig. 8. Frequency response of the magnetometer, mea- 
sured by recording the response to an oscillating mag- 
netic field generated by a test coil. The dots show mea- 
sured points recorded in free running mode under condi- 
tions optimized for maximal magnetometric resolution, 
(a) Calculated first-order low-pass filter corresponding 
to a spin polarization lifetime of 1.67 ms. (b) Fitted 
frequency response taking into account (a) and the 4**^ 
order low-pass filter of the lock-in amplifier (time con- 
stant = 30 lis), (c) Measured frequency response in the 
phase-stabilized mode. 



for which both approaches yield the same magnetometric 
sensitivity. For bandwidths larger than /q, a phase mea- 
surement is more advantageous whereas for bandwidths 
smaller than /o a frequency measurement gives the higher 
sensitivity. 

B. Bandwidth 

In addition to the sensitivity, the bandwidth, i.e., the 
speed with which the magnetometer signal follows mag- 
netic field changes, is an important feature of a magne- 
tometer. The steady-state solutions of the Bloch equa- 
tions, Pqu and Pip [Eqs. Q and |0J], follow small field 
changes at a characteristic rate r2, corresponding to a 



delay time taus 



Since the steady state is only 



reached exponentially, the frequency response is that of 
a first order low-pass filter [see Fig. Oa)] with a (-3 dB) 
cut-off frequency fc given by 



JC == -X = T^ = Al/2 , 



and hence a bandwidth of 
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where Ai/2 is the half width of the phase signal measured 
in Hz. 

To achieve maximum sensitivity, atomic magnetome- 
ters aim at a maximum Tg , at the cost of a reduced band- 
width of typically a few tenths of Hz. A large bandwidth 
can be obtained by increasing the light power since that 
leads to shorter ts and therefore to higher bandwidth. 



Larger light powers also increase the S/N ratio but the 
effect can be overcompensated by magnetic resonance 
broadening, resulting in a degradation of the magneto- 
metric resolution. 

Using feedback to stabilize the magnetic resonance 
conditions is another way to increase the bandwidth. 
Figure |SJc) shows the frequency response of the 0PM 
in both the free-running (without feedback) mode and 
in the phase-stabilized mode where the phase signal is 
used to stabilize Lo-d to the Larmor frequency lo]^. For 
large loop gain the bandwidth is mainly limited by loop 
delays. 

A third method to achieve large bandwidths is the so- 
called self-oscillating mode. In this mode the oscillating 
signal measured by the photodiode is not demodulated 
but rather phase-shifted and fed back to the rf-coils. For 
a 90° phase shift the system then oscillates at the Larmor 
frequency. In order to measure the magnetic field, the 
frequency of this oscillation has to be measured. Mag- 
netic field changes then show up — at least theoretically^ 
— as instantaneous frequency changes. In practice, reac- 
tion times smaller than a single Larmor period have been 
observedl^. 

Of the three modes outlined above, the latter two both 
rely on frequency measurements. The self-oscillating 
magnetometer provides a frequency that has to be mea- 
sured. The phase-stabilized magnetometer measures the 
frequency via a reference frequency locked to the Lar- 
mor frequency. As a consequence, both methods suffer 
from the reduced magnetometric resolution predicted by 
Eq. H18I) . Therefore, we have concentrated on the free- 
running mode of operation for which the magnetometric 
resolution is given by Eq. 119|) and the bandwidth by 
Eq. IIH). 

Thanks to the rather high light power required for 
optimal magnetometric resolution at higher cell tem- 
peratures, the cut-off frequency of the magnetometer 
was 95 Hz. The bandwidth of the device under these 
conditions can be extracted from the transfer function 
(Fig.|HJb)) and is about 140 Hz. Because of the time con- 
stant of the lock-in amplifier, the measured bandwidth is 
10 Hzsmaller than the 7r/2 x 95 Hz one would expect for 
a first order low-pass filter [Eq. 1221) ]. 

C. Experimental lineshapes 

Figure 1^1 shows a Nyquist plot with experimental data 
and a model simultaneously fit to the in-phase and 
quadrature components of the data. The data show a cer- 
tain asymmetry that can not be reproduced by the model. 
The Nyquist plots for different magnetic field distribu- 
tions (Fig. [SJ suggest that the asymmetry is caused by 
inhomogeneous magnetic fields. Unfortunately the mod- 
els discussed in Sec. 12 Bl do not fit the data correctly, im- 
plying that higher-order gradients cause the deformation 
of the measured lineshape. The fact that the asymmetry 
is more pronounced for high rf amplitudes indicates that 
inhomogeneous rf-fields — causing the different parts of 
the ensemble to contribute with different widths — have 
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Fig. 9. The dots represent a Nyquist plot of the magnetic 
resonance measured under optimized conditions. The 
sohd hne is a fit of Eq. (|15|l with added offset and phase 
rotation to the measured data. The fit model assumes a 
constant magnetic field distribution. The offset is indi- 
cated by the dot close to the origin. The short diameter 
of the ellipse is drawn in order to illustrate the phase 
rotation of 2.4°. 



to be considered. Unfortunately, models for such inhomo- 
geneities do not lead to analytic line shapes. An empirical 
model which assumes the measured resonance consists of 
a sum of several resonances, each at a different position 
and with a different width, can be fit to the data. The 
data can be fit perfectly if the number of subresonances 
is high enough. However, such fits have a slow conver- 
gence and do not provide the needed information about 
the width and amplitude of the resonance in single fit 
parameters. For practical reasons (during the optimiza- 
tion more than 2000 spectra were fit) we decided to use 
the constant magnetic field distribution model for fitting 
data similar to the ones in Fig. |51 

Magnetic field inhomogencities have much less influ- 
ence on the shape of the phase signal resulting in more re- 
liable values for r2 . The phase signal represents the speed 
with which the resonance evolves through the Nyquist 
plot. Using both the phase signal and the Nyquist plot, 
the in-phase and quadrature components of the resonance 
were reconstructed, however, the frequency scaling were 
given by the phase signal only. 

D. Optimization measurements 

For the optimization of the NEM given by Eq. ifT^I) the 
S/N ratio of the lock-in input signal and the linewidth r2 
have to be measured. Figure [T] shows a frequency spec- 
trum recorded at the input of the lock-in amplifier using a 
FFT spectrum analyzer. The frequency Wrf was tuned to 
the center of the magnetic resonance so that the modula- 
tion of the photocurrent was at its maximum amplitude. 
The power spectrum shows a narrow peak at Wrf sur- 
rounded by noise peaks that characterize the magnetic 
field noise. Monochromatic magnetic field fluctuations, 
e.g., line frequency interference, modulate the phase of 



the measured sine wave and show up in the power spec- 
trum as sidebands. The low frequency flicker noise of the 
magnetic field thus generates a continuum of sidebands 
that sum up to the background structure surrounding 
the peak in Fig. [T] The estimation of the intrinsic sen- 
sitivity is based on the assumption that those sidebands 
would disappear in a perfectly constant magnetic field. 
The amplitude noise of the signal is mainly due to the 
electron shot noise in the photodiode, which generates 
a white noise spectrum. For frequencies which are more 
than 1 kHz away from the resonance, the noise level drops 
to the white noise fioor. The electron shot noise is the 
fundamental noise level that can not be avoided. The 
noise spectral density pi can be calculated from the DC 
current /dc flowing through the photodiode: 



PI 



V2e/DC 



(23) 



At room-temperature the measured rms noise in the 
spectrum was 5% to 20% above the shot-noise level, de- 
pending on induced noise on the photocurrent and the 
laser frequency stabilization that could cause excess noise 
in the light intensity. The rms noise rose rapidly for 
higher temperatures because of the increasing leakage 
current in the photodiode. Unfortunately, in the ex- 
perimental setup the photodiodes were in good thermal 
contact with the Cs cell and, given that the optimal op- 
erating temperature of the Cs cells turned out to be in 
the range of 50° C to 60° C, the photodiode produced 
an excess noise larger than the shot noise of the pho- 
tocurrent. Figure shows a spectrum recorded under 
conditions optimized for maximal magnetometric resolu- 
tion. At 53° C the measured rms noise was higher than 
the shot noise by a factor of 1.55. However, this limita- 
tion can be overcome easily since the photodiodes do not 
need to be close to the Cs cell and thus can be operated 
at room-temperature. In order to avoid the problem of 
drifting values of p during the optimization of pB, the 
theoretical shot noise level pi was used for p in Eq. El 
instead of the measured noise. 

The amplitude A of the signal can be extracted from 
the FFT-spectrum by integrating the spectrum over 
three points (±1 Hz) around the center frequency. The 
procedure was needed since the Hanning window used 
by the spectrum analyzer to reconstruct the spectrum 
causes a slight broadening of the central peak. The val- 
ues calculated in that way are in good agreement with 
those measured by the lock-in amplifier. 

The third parameter needed to calculate the intrinsic 
sensitivity is the half-width T2 of the magnetic resonance. 
The value was extracted from a magnetic resonance spec- 
trum recorded by the lock-in amplifier during a frequency 
sweep of the applied oscillating magnetic field. As dis- 
cussed in Section 14 CI a constant-gradient model was fit 
to the data in order to extract r2. 

For optimizing in a three-dimensional parameter space, 
the time for one measurement had to be kept as short as 
possible. When the lock-in amplifier signal was used as 
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Fig. 10. Dependence of the magnetic resonance 

linewidth on the rf amplitude C/rf. The points are ex- 
tracted from measured magnetic resonance spectra by 
least squares fitting of model Eq. (fT^ . The phase sig- 
nal (a) has a constant linewidth whereas the common 
widths of the in-phase and quadrature signals (b) in- 
crease rapidly with rf amplitude. The solid line repre- 
sents a model fitted to the data that assumed an ad- 
ditional broadening caused by inhomogeneous magnetic 
fields. 



a measure for A [see Eq. (fT7|) ] and the noise was calcu- 
lated from the DC current it was not necessary to record 
a FFT spectrum for each set of parameters of the op- 
timization procedure. In that way the time for a single 
NEM measurement was reduced to the 20 s sweep time of 
LUrf plus the time needed to measure the DC current and 
the temperature of the cell. The measurement was con- 
trolled by a PC running dedicated software for recording 
and fitting the magnetic resonance signals. The ampli- 
tude of the rf field, B^-f, was changed automatically by 
the software, resulting in series of typically ten NEMs as 
a function of B,-{. A typical optimization run was made 
by recording many such series while the system slowly 
heated up. Repeating those runs for different light pow- 
ers finally resulted in data for the whole parameter space. 

5. Results 

A. Dependence on rf amplitude 

The first study made with the magnetometer examined 
the dependence of the magnetic resonance on the rf am- 
plitude i?rf measured a series of spectra recorded at room 
temperature. Figure El shows the dependence of the 
magnetic resonance signal width on the rf amplitude 
measured by the coil voltage C/,.f . The width of the phase 
signal (see Fig. IIOII was fit with a constant, whereas the 
common width of the in-phase and quadrature compo- 
nents were given by Eq. Q. To fit the widths correctly, 
a constant width had to be added to Eq. (0. The ad- 
ditional constant width can be interpreted as a residual 
broadening caused by magnetic field inhomogeneities of 
higher order than the one considered in the line fitting 




Fig. 11. Amplitude of the in-phase (ai) and quadra- 
ture (02) signals as a function of rf amplitude C/rf. The 
points represent values extracted from measured mag- 
netic resonance spectra. The solid lines show a model fit 
to the data points (see text) . The quadrature amplitude 
a2 is equal to the amplitude of the incoming sine wave 
on resonance {S = 0). 



model. The Nyquist plot (see Fig. O shows that higher 
order gradients are present and the excellent agreement 
in Fig. [TUl suggests that they can be treated as an addi- 
tional broadening. 

Figure 1111 shows the amplitudes of the in-phase and 
quadrature magnetic resonance signals. The amplitudes 
where extracted from the same spectra used for Fig. 1101 
The fit model used to explain the amplitudes (solid lines 
in Fig. lll|) was based on Eqs. H1U|) and (|ll|l with a back- 
ground proportional to B,-{. The origin of the background 
was an inductive pick up of the Bi field by the photocur- 
rent loop which caused an additional phase-shifted sine 
wave to be superposed on the photocurrent. As discussed 
in the theory part (see Fig. 2)) that lead to an offset in the 
measured amplitudes of the magnetic resonance signal. 

The NEM as a function of rf amplitude is inversely 
proportional to the quadrature amplitude (02 in Fig. lll|l . 
since the linewidth of the phase signal and the shot noise 
do not change with rf amplitude. The optimal rf ampli- 
tude was determined from the data shown in Fig. II II and 
corresponds to S* = 1. 

B. Dependence on temperature and light power 

As described in section 14 PI the dependence of the NEM 
on the temperature was recorded while the system was 
slowly heated. The rf amplitude was automatically 
scanned so that for every temperature the optimal rf am- 
plidudc could be determined. Figure [T^ shows a contour 
plot of the NEM as a function of temperature and light 
power. If the light power is increased, the temperature 
(and hence Cs atom density) has also to be increased 
to maintain optimal resonance conditions. Figure ^Jb) 
shows the power transmitted through the cell relative to 
the incident light power. A relative transmission of 0.37 
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For light powers below 15 fiW and the corresponding 
temperatures the sensitivity rapidly degrades. The loss 
in sensitivity is less pronounced if the power and tem- 
perature are chosen above the optimum. Values for i?gN 
of up to 500 000 (114 dB) were measured at a resolution 
bandwidth of 1 Hz. 

The optimal magnetic field resolution of our magne- 
tometer is reached at a light power of 15 fiW and a tem- 
perature of 53° C. With that set of parameters, the us- 
able bandwidth of the magnetometer was determined by 
a cut-off frequency of about 80 Hz (see Fig. Iisr a)). In 
order to meet the required 100 Hz bandwidth a slightly 
larger light power can be used. All characterizing mea- 
surements (cf. Figs. 12 |S1 El and [TJ were therefore per- 
formed with a light power of 20 /iW at 54° C. 



Fig. 12. Contour plot of the magnetometric resolution 
(NEM) as a function of temperature and light power. 
The map is calculated by numeric interpolation from 
seven optimization runs (indicated by vertical lines). The 
labels at the contours mark the NEM in fT/Hz^/^. The 
points of minimal NEM for each optimization run are 
indicated by points. The connecting line is a cut along 
which the data of Fig. E| are obtained. Including the 
variation of the rf amplitude 970 parameter sets were 
recorded and analyzed to produce the map. 
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Fig. 13. Magnetic resonance parameters as a function 
of incident light power: (a) shows the width of the phase 
signal that determines the cut-off frequency of the mag- 
netometer bandwidth; (b) shows the DC transmission 
through the Cs cell relative to the incident light power; 
(c) is the signal/noise ratio of the lock-in input signal 
with respect to calculated shot noise; (d) shows the NEM 
at the points indicated in Fig. El 



corresponds to an absorption length which matches the 
cell length. Taking into account losses at the windows, 
a density corresponding to 1.4 absorption lengths was 
found to be optimal. 

For each light power the optimal temperature is indi- 
cated by a dot in Fig. [HI Plotting the NEM along the 
optimum temperature power curve, i.e., the curve con- 
necting the dots, results in the plot shown in Fig. 113r d). 



6. Conclusion 

Optimizing the performance of the magnetometer has 
led to a set of parameters for which the device offers 
a large sensitivity and a large bandwidth. Both require- 
ments can be met at the same time because of rather 
large linewidths that turned out to be optimal. Under 
these conditions the high magnetometric sensitivity re- 
lies on the achieved very high signal/noise ratios. The 
system has the potential to operate at a i?sN of 500000 
(Fig. O and we hope to be able to demonstrate this once 
the photodiodes can be removed from the heated Cs cell. 
However, even using the measured i?sN of 320000, the 
intrinsic sensitivity of 100 fT/Hz^ " is good enough for 
less demanding cardiomagnetic measurements. 

The magnetometer bandwidth of 140 Hz in the free- 
running phase-detecting mode (Fig. jSl is high enough 
for cardiac measurements. The phase-detecting mode 
avoids the fundamental limitations associated with fre- 
quency measurements using short integration times. 

An important open experimental question is whether 
the predicted intrinsic sensitivity can be reached using 
several of the present OPMs in a higher order gradiome- 
ter geometry. With gradiometric SQUID sensors it is 
possible to achieve NEM value on the order of 20fT/Hz^ " 
in unshielded environmentsi&. In future we plan to 
use cells with spin-preserving wall coatings rather than 
buffer-gas cells as sensing elements. Coated cells have 
the advantage that the atoms traverse the volume many 
times during the spin coherence lifetime, therefore aver- 
aging out field inhomogeneities. We are therefore confi- 
dent that the present limit from field gradients can be 
overcome and that optical magnetometers can reach an 
operation mode limited by their intrinsic sensitivity. 

7. Appendix A: The Cramer— Rao bound for fre- 
quency measuring magnetometers 

For the measurement of the frequency w of a sine wave 
with a rms amplitude A sampled at TV ^ 1 points sepa- 
rated by time intervals Tg the Cramer-Rao lower bound 
for the variance V^ of w in the presence of white Gaussian 
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amplitude noise of variance cr^ is given hy^: 



K, = 



12a' 



A'^TljN' 



(24) 



where Tm — NTs is the total time interval for one 
frequency determination. The bandwidth on the input 
side of the lock-in amplifier is therefore Fhw = l/^Ts = 
N/2Tm, that at the output is /bw = l/2TAf. With the 
definition of the signal-to- noise ratio, Eq. (|17|l . \{j can be 
expressed independently of the number of samples: 



K.= 



12p?Fi 



bw 



6 



A'Ti,N 



d2 7^3 ■ 



(25) 



Ideal measuring processes are limited by that condition 
only. Frequency measurements by a FFT with peak in- 
terpolation is a Cramer-Rao bound limited measuring 
process-'*. 

From that bound a lower limit for the performance of a 
frequency measuring magnetometer can be derived. The 
so-called self-oscillating AI^ magnetometes^ is of this type 
since it supplies an oscillating signal with a frequency 
proportional to the magnetic field. With Eq. (QJ it follows 
that the root spectral density of the measurement noise 
Pb is given by: 



PB = 



/bw 7F-RsN 



(26) 



8. Appendix B: The Cramer Rao bound for 
phase measuring magnetometers 

The Cramer-Rao lower bound for the measurement of 
the phase of a signal with known frequency is given hy^: 



K, 



A^N 



(27) 



An example of a measurement process limited only by 
that condition is the lock-in phase detection where the 
phase is calculated from the in-phase and quadrature out- 
puts of the lock-in amplifier [see Eq. Q]. In order to 
calculate the variance V^ of the phase measurement we 
assume a white amplitude noise spectrum with a power 
spectral density p^: 



K. 



A^2T, 



M 



(28) 



Using this expression and the definition of i?sN [Eq. lfT7|) ]. 
Eq. H28f) can be written as 



= V^ = 



P /bw 

A^ 



/i 



bw 



p2 



(29) 



From the measured phase (p, the detuning 5 — w^f — 
lul can be derived. For ip ^ 1, Eq. (TJ leads to ^ « 
r2(p- Using Eq. Q, the detuning can be expressed as a 



magnetic field diff'erence AB = 5/j which leads, together 
with Eq. (|29|l . to the magnetic field resolution as- 
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If 



CTipF 



rv/bw 



(30) 



7f 7F-RsNv7b^ 
The root spectral density of the noise in the As mea- 
surement, pB = <7b/ f\^!„ , is thus given by: 



PB = 



(31) 
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